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Extending the U(l) slave-rotor representation[l] of the repulsive Hubbard model, we propose an 
SU(2) slave-rotor decomposition for the attractive Hubbard model, where the SU(2) slave-rotor vari- 
ables represent order parameter fluctuations associated with superconductivity and charge density 
wave. This decomposition method allows us to modify the standard Hartree-Fock mean field theory 
by incorporating order parameter fluctuations on an equal footing. Deriving an effective SU(2) 
slave-rotor action from the attractive Hubbard model, and analyzing it at the mean field level, we 
demonstrate a second order phase transition driven by softening of the slave-rotor variables. 
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Recently, Florens and Georges proposed a slave-rotor 
representation decomposing bare electrons Ca into col- 
lective charge excitations e~'^ and renormalized elec- 
trons /cr, i.e., Ccr = e~'^^fa-[l] Applying this represen- 
tation to the repulsive Hubbard model, they explained 
the Mott-Hubbard transition from a spin liquid Mott 
insulator to a Fermi liquid metal in the square lattice. 
Although the basic scheme of the rotor representation 
is quite appealing, [1-3] there are several unsatisfactory 
points in this approach. First of all, effects of spin fluc- 
tuations are not well described. This is the reason why 
the previous slave-rotor theories considered only param- 
agnetic phases. [1-3] Even if only charge fluctuations are 
taken into account, the slave-rotor representation is not 
complete in the sense that the SU(2) pseudospin sym- 
metry in the Hubbard model[4] is not reflected in the 
U(l) slave-rotor representation. Recently, we extended 
the U(l) rotor formulation into an SU(2) one for SU(2) 
charge fluctuations in the repulsive Hubbard model. [5] In 
this study we found an anomalous metallic phase with a 
pseudogap. 

In this paper we apply the SU(2) slave-rotor 
representation[5] to the attractive Hubbard model. One 
main difference from the study of the repulsive Hubbard 
model is the presence of nonzero order parameters asso- 
ciated with superconductivity (SC) and charge density 
wave (CDW). An important task is to develop how to 
incorporate order parameter fluctuations in the conven- 
tional mean field theory such as the BCS scheme. A 
standard weak coupling procedure is to integrate out 
electron excitations and expand the resulting logarithmic 
action for order parameter fluctuations around the mean 
field ground state, called the Landau-Ginzburg- Wilson 
(LGW) expansion. Although this approach is system- 
atic and firm-based, the procedure to integrate out gap- 
less electrons near the Fermi surface gives rise to several 
uncertainties in the LGW effective theory, [6] thus it is 
necessary to treat both electron excitations and order 
parameter fluctuations on an equal footing. 

In this paper we derive an effective theory from the at- 
tractive Hubbard model, imposing order parameter fluc- 
tuations at the mean field level. Our mean field action 
consists of two parts. One is a fermion sector correspond- 



ing to a modified Hartree-Fock theory, and the other a 
boson part reflecting order parameter fluctuations, thus 
allowing us to take into account both electrons and order 
parameter excitations on an equal footing. It turns out 
that SU(2) rotor variables represent order parameter fluc- 
tuations, and their presence in the effective action admits 
us to analyze effects of their fluctuations in the mean field 
approximation. We discuss how this mean field scheme 
modifies the conventional Hartree-Fock theory. 
We consider the attractive Hubbard Hamiltonian 

H ^~tY^ {c\^cj„ + H.c.) - ^-Y^c\^c^^c\^c^lll) 

(i])a i 

where t is a hopping integral of electrons, and g a coupling 
constant of effective attractions. The interaction term 
can be decomposed into pairing and density channels in 
the following way 

_^ t t _ _5 t t _ t t 

Using the identity 

n^^nil = i(ni| + n,,i - if + ]^{n,^ + n,i - 1) 

with riicr — c\^Cia, and performing the standard 
Hubbard-Stratonovich (HS) transformation for each in- 
teraction channel, we obtain 

Z = J Dc,,D^,D^,e- fo''^'', 

icr {ij)a- 
i i 
i rr " i 

Here (jji is an SC order parameter associated with an ef- 
fective pairing potential, and ipi is a CDW order parame- 
ter involved with an effective density potential. We note 
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that the chemical potential /x = fj,b + g/2 differs from its 
bare value /ifc. 

One can represent Eq. (2) in terms of a Nambu spinor 
for convenience in describing superconductivity. [7] Using 



the Nambu spinor tpi = ( 1 ) '^^^ obtain 



Z = J DilJiD<j)uDcl>2iD<j,sie-^odrL^ 



(ij) 



i 



(3) 



where (pu and (j)2i are the real and imaginary parts of the 
pairing potential 0i, i.e., (f>i = (pu — i4>2i, and (fii is re- 
placed with for a unified notation of SU(2) symmetry. 
Introducing a pseudospin variable Vli = ((/iii, <^2i, f/'si), 
one can express Eq. (3) in the following compact form 



(ij) 



(4) 



Since the above effective Lagrangian is quadratic in 
electron excitations, one can formally integrate out the ■i/'j 
fields to obtain an effective Lagrangian of the pseudospin 
order parameter by expanding the resulting logarithmic 
term for the fij fields. It should be noted that the expan- 
sion parameter is g / D. where D is an electron bandwidth, 
thus this expansion can be justified in g/D « 1. Al- 
though the weak coupling condition is satisfied, there are 
still several unsatisfactory points in this order parameter 
action. It is difficult to justify the LGW expansion in 
the presence of gapless electrons because they can cause 
nonlocal interactions between order parameters, making 
it unreliable a conventional treatment in a local effec- 
tive action. [6] In this respect we do not integrate out 
the fermion excitations in deriving an effective action. 
Instead, we treat both electrons and order parameter ex- 
citations on an equal footing, as mentioned before. 

For the equal treatment of electrons and order param- 
eters, we apply a strong coupling approach to this prob- 
lem, meaning to solve the interaction term first. Remem- 
ber that the weak coupling approach is to solve the ki- 
netic energy term first and treat the interaction term per- 
turbatively based on the non-interacting fermion ensem- 
ble. Using the identity called the CP^ representation [8] 



. f=miUiT3UL 



Ui = 



, (5) 



where Zi^ in the SU(2) matrix field Ut is a boson field 
with pseudospin a to satisfy the unimodular constraint 
Scr N»<^P — 1' performing the gauge transformation 



(6) 



one can solve the coupling term from — miiljl{(li-T)^i 

to — Ylii "^i??!' T'3'?i! where vrii is an amplitude of the pseu- 
dospin order parameter. We call Zi^ and 'qia a bosonic 
spinon and a fermionic chargon, respectively. 

In this strong coupling approach we find an interest- 
ing physics that order parameter fluctuations $1; carry- 
ing a pseudospin quantum number 1 fractionalize into 
bosonic spinous Zi^ with pseudospin 1/2 in order to 
screen out the pseudospin of an electron due to strong 
interactions. The components of the rji field are given 



t 



-ZilCi^ + Zifcl^ 

say this fractionalization is that bare electrons ipia- frac- 
tionalize into bosonic spinous Ui^ra' and fermionic char- 
gons r]icr, i.e., tjjic, = Uiccr'Tjicr' owing to strong interac- 
tions. 

Inserting Eqs.(5) and (6) into Eq. (4), we obtain 
Z = J Dr]iDUi5{U}Ui-l)e-^o'irL^ 



Another way to 



-tY^ivlulrsUjrij+H.c.) 



(ij) 



(7) 



Note that the integration measure / DilJiDCli in Eq. (4) 

is changed into / Dri,DU^5{ulu^ - 1) in Eq. (7). The 
number of integration variables is 4 in both cases. One 
important point in this expression is that we impute in- 
teractions between electrons and order parameter fiuctu- 
ations to couplings between chargons (renormalized elec- 
trons) and spinous (fractionalized pseudospins) in the ki- 
netic energy. A standard way to treat this nontrivial 
kinetic energy term is to integrate out the chargon fields. 



Z = j DUi5{U}Ui - 1) exp [- dr^ ^ t^KTs) 
+trln(M - m^T3 - iiU^TsUi + UjOrUi - tUjraU, 



(8) 

where tr in the logarithm means sum over time, space, 
spin and matrix elements. Expanding the logarithmic 

term for the bosonic spinous Ui^a' (zia-), one obtains an 
effective action of the spinous. One important difference 
from the weak coupling approach is that the expansion 
parameter is D/g instead of g/D. Unfortimately, this 
conventional strong coupling approach has an important 
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defect. Metallic physics (information of a Fermi surface) 
of electrons is not introduced in this expression. Actu- 
ally, expanding the logarithmic term in the expansion 
parameter D/g, the resulting effective action is known to 
be the 0(3) nonlinear u model appropriate to an insu- 
lating magnet, describing competition between SC and 
CDW.[9] 

Instead of integrating out fermions, we decouple the 

"interacting" kinetic energy into the conventional "non- 
interacting" one via the HS transformation 



(9) 



where Eij and are HS matrix fields associated with 
hopping parameters of the rji fermions and Ui bosons, 
respectively. They arc sclf-consistently determined from 
the saddle point equations in the mean field approxima- 
tion 



E}f = ivlv.s), = {Ulf,Tsp^U,^s). (10) 

We make an ansatz for the hopping matrices as 

E^' = Etsc.j, F,f = FTsc.j, (11) 

where E and F are amplitudes of the hopping parame- 
ters. The reason why we introduce the T3 matrix is that 
the fermion sector should recover the original electron 
Lagrangian Eq. (4) as the slave-rotor representation [3] 
does. 

We also perform a mean field decomposition in the 
coupling term of the time part as 

with the mean field ansatz of 



(13) 



This ansatz is consistent with Eq. (11). The chemical 
potential term is decoupled in the mean field level as 

= - I^VialiaSViS - l^hasU}^i^T30^UijS + jJ-kaSliaS (14) 



with 



QiaS 



(15) 



Inserting Eqs. (9), (11), (12), (13), (14), and (15) into 
Eq. (7), we find an effective Lagrangian for the mean 
field analysis of the attractive Hubbard model 



drL 



Z = J Dr]iDUiS{UjUi - l)e-^o 
L = Lq + + Ljj, 

Lo = 4t 5^ £F - 2 ^ + 2/x 5^ gJi + ^ m?, 

(ij) i i ^ i 

i 

-tFj2in}T3Vj + H.c), 

(ij) 

Lu = Y^^MUldrUiTz) - tEY,tv{UlTzUjT3 + H.C.) 



(ij) 



(16) 



In this effective Lagrangian qi and TOj are not all 

independent. One can easily sec 2li = m,i/g. Introducing 
—hri = hi — jiqi in Eq. (16), we obtain 



Z = J D7^,DU,d{UjU, - l)e-^o<i^L, 

L = Lq + Ltj + Ljj, 

Lo = Uj2EF+^^hrimi + j^Yl 



(ij) 



= ^'ni{dA - ruiTs - hriT3)r]i - ^-^^(TytTsT?^ + H.C.), 



(ij) 



Lu = ^Yl mitr{U}drUiT3) -tEY^ tr{U}T3UjT3 + H.c.) 

9 ■ 



{ij) 



-^Y,mMUhsUir3). 



(17) 



To explore the consequences of order parameter fluc- 
tuations in the effective rotor action Eq. (17) at the 

saddle point level, wc express Eq. (17) in terms of the 
spinous Zia- with the mean field ansatz rrii = (— l)*m and 
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Z = J DrjiDzi^e-^o'i^L, 

L = Lq + Lri + Lz, 



Lr, = J2vl{drl - (-l)'mr3 - (-l)'/i.r3)r?i 

i 

-tFj2ivlT3rjj+H.c.), 

{ij) 

i<7 {ij)cr 

_^^j2{-iy<rzU. + aE(E i^-i' - 1)' (18) 

where A is a Lagrange multiplier field to impose the uni- 
modular constraint, replaced with its mean field value. 
We also replaced 2E with E in the above. 

The SU(2) slave-rotor action consists of two parts. One 
is the fermion sector L^, and the other the boson part 
Lz- The fermion action coincides with a conventional 
mean field theory, the Hartree-Fock theory except for 
the renormalization of the bandwidth t tF , \i we re- 
gard m -f /ir as an effective magnetic field. At zero tem- 
perature the fermion sector thus always remains in the 
" magnetic" phase with a gap to quasiparticlc excitations 
due to Fermi-nesting if half filling in the square lattice is 
considered. However, it is important to see that quan- 
tum fluctuations of the Zi^ bosons reduce the staggered 
pseudo-magnetization m, compared to the Hartree-Fock 
magnetization. The boson action, on the other hand, 
can be considered as corrections to the fermion mean 
field theory due to order parameter fluctuations, not cap- 
tured in the LGW framework. It is basically the same as 
the CP^ Lagrangian of the 0(3) nonlinear g model, but 
modified by the Berry phase term y J^iai^^T-^lu^-rZia- 
The presence of tr = ± in the spinon Lagrangian L^, 
resulting from the T3 matrix, also leads to modification 
of the CP^ Lagrangian. If the T3 matrix is not utilized 
in Eq. (11), the hopping term in vanishes, thus Eq. 
(4) cannot be recovered from Eq. (18) by following its 
reverse procedure. 

In analyzing the effective SU(2) rotor Lagrangian Eq. 
(18), we confine our attention to half-filling (/x = 0) in 
the square lattice for simplicity, and discuss effects of hole 
doping later. Let us begin with the fermion sector. The 
fermion Lagrangian is well known because its structure is 
nothing but the BCS theory. The mean-field conditions 
for m and £ at T = read 



2g{m + hr) 



Ft 



V2 



(19) 



Here e]! = — 2t(cos fc^. -l-cos ky) is the bare band in the ab- 
sence of effective exchange splitting introduced by non- 



zero m + hr, and E^ = ^ (Ee^)'^ + {m + h^Y is the char- 
gon energy with a gap set by m 4- ^r- D = At is a half 
of the bandwidth, and the fc-sum in both equations is 
over the reduced Brillouin zone. One important differ- 
ence from the conventional Hartree-Fock scheme is the 
presence of hr, modifying the first equation. 

The spinon Lagrangian is given in the energy- 
momentum space 

= E ^Iku + Eael) Zaku 

aki/ 
I 

+ E ^Ik+Qu - Eael) Zak+Qu 

akf 

I 

+ — Z^l'^y^ak+Qv^'^ku + Z^k^Zak+Qu) 



crkv 



'Y'y{zl^j^Q^Zakv + zlj.^Zak+Qv), (20) 
^ (7kv 

where e| = —2t{coskx + cos ky) is the bare spinon dis- 
persion with Q = (7r,7r). The bosonic fc-sum is also over 
the reduced Brillouin zone. The boson part can be diag- 
onalized using a pair of operators {'yc7ikv,1i72kv) related 

to {Zaku,Zak+Qv) by 

Zaku = --^(cosh6l<^fc - sinh^<^fc)(7<^ifcj. +7^2fei/), 

Zak-i-Qi^ = -^{coshe^k + Smh0^k){lalkiy - 7<72ki^pl) 

After taking cosh26'^A; = A/£|, sinh2^o-fe = Eael/E^, 
and El = ^J'>?^^^{Ee^, one gets 



= *7 E ^'Clllkullku - l]y2kul<T-2ku) 

/xm^ I t ^ 

- — 2^ Cr(7alfci.7lfci^ - 7a2fei.7^2fc<.) 
akv 

t 

+ E ^killlki.l'r-i^ku + 7l2fe^7a2fe^). (22) 
kv 

The boson spectrum E^. is gapped if A — ED > while 
A — ED — leads to the condensation of Zia ■ 

The mean field equations for A, F, and hr are obtained 
to be 



-e4 df=y. 



Eg el; 



z2 



5A ^E, 



k ^^'k 



and 



13^ Ef + u^7v?/g^■ 



{22.) 



(24) 



An important point in Eq. (24) is that if we do not in- 
troduce an energy cutoff in the frequency integral, the 
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mean field equation (24) diverges after doing tlic fre- 
quency integral. However, one need not be surprised at 
this divergence because it also happens in the mean field 
equations for E and F unless we introduce a momen- 
tum cutoff, here the bandwidth D. In this respect it is 
necessary to introduce an energy cutoff. It is natural to 
take the energy cutoff as m because the presence of m al- 
lows this decomposition. When evaluating the frequency 
integral, we first divide the integral into two parts, diver- 
gent and divergent-free parts. We calculate the divergent 
part within the energy cutoff, but the divergent-free part 
without the energy cutoff, i.e., performing the Matsub- 
ara summation in the divergent-free integral. As a result, 
Eq. (24) reads 



2ff 9^ V- 



(25) 



To get an idea on the analytical structure of the set of 
self-consistent equations obtained above, we first rewrite 
Eqs. (19), (23) and (25) as the integration over the en- 
ergy with a certain density of states Die), and approx- 
imate it with a constant value, D{e) = 1/{2D), thus 

J2'k — (1/4^) J^jj de. The mean-field equations are then 

given by 



m + hr = 



FD 



tFD 



E 



sinhf^^ i-T— I 

smh[2FD/g] - 2FD/g 



2mED 



5A 



Asm\i[FD/g] 
= sm-^{ED/X), 



F = 



AmED/g\ - sm[AmED/g\] 
^{m/g)sm^[2mED/g\\ ' 



IT 



AED 



(26) 



The Bosc condensation occurs at Ac = EcD, giving rise 
to nic/g — 7r/4 K. 0.8. The critical {g/D)c is determined 
from the following equation 



1 jP/g) 

4 



[sinh(Z?/5), - (D/g), 



1 



sinh^[(l/2)(i?/.g)e] 

{D/g)c 



' sinh 



(^/2){D/g), 



Tr/i+2/v- 



(J)/g)c|Binh(J)/g)c-(D7g)^ 
ih^[(l/2)(D/g)cl 



_ TT 2 

4 tt' 



(27) 



Solving this equation numerically, we find {g/D)c ~ 0.3. 
This means that the Bose condensation of spinous ap- 
pears in g/D < {g/D)c while they become gapped in 
g/D>{g/D),. 

When the bosonic spinous are condensed, this phase is 
identified with an SC state or a CDW phase. We cannot 



determine which phase arises because of the SU(2) pscu- 
dospin symmetry in the Hubbard model at half filling. If 
easy plane anisotropy is introduced, the resulting state 
would be an SC phase. When the spinous are gapped, 
the resulting phase is an insulating state owing to a gap 
in the chargon spectrum at half filling. Remember that 
chargon excitations are always gapped at half filling due 
to Fermi-nesting. Thus, the SU(2) slave-rotor mean field 
theory shows a second order phase transition from an SC 
state to an insulating phase at half filling, varying the 
strength of local attractions. One may interpret the in- 
sulating phase as a preformed-pair state due to strong 
phase fluctuations of preformed pairs. 

So far, we performed a saddle-point analysis and ob- 
tained a mean-field picture, showing a second order phase 
transition for the spinon field Zi^ associated with order 
parameter fluctuations. It is natural to ask the stability 
of the mean-field picture against gauge fluctuations aij 
that appear in the phase fluctuations of the hopping pa- 
rameters, Eij = Ee^^^^^^Ts and = Fe^'^'^'^^Ta, where 
E and F are the mean field values obtained before. It 
should be noted that the U(l) pseudospin-gauge field 
is compact, thus allowing instanton excitations [10]. From 
the seminal work of Fradkin and Shenker[ll] we know 
that there can be no phase transition between the Higgs 
and confinement phases. The order parameter discrimi- 
nating the Higgs phase from the confinement one has not 
been known yet. In this respect only a crossover behavior 
is expected. In the present problem the phase-coherent 
state corresponds to the Higgs phase while the phase- 
incoherent state coincides with the confinement phase. 
Applying Fradkin and Shenker's result to the present 
problem, we conclude that the second order phase tran- 
sition turns into a crossover between the coherent and 
incoherent phases. The pseudospin order parameter, be- 
ing a gauge-invariant quantity, remains unaffected by the 
gauge fiuctuation. 

One cautious person may ask what the crossover means 
physically. Since a systematic method for evaluating the 
electron spectral function in the confinement phase is not 
known, it is difficult to say any physical statements for 
the spectral function in fact. However, recalling the re- 
cent study that considerable portions of coherent quasi- 
particle spectral weights are transferred to incoherent 
backgrounds as the local interactions are increased, [12] 
the crossover in this paper can be understood by trans- 
ferring the coherent spectral weights of quasiparticles to 
the incoherent backgrounds. Actually, our mean field 
analysis coincides with this picture. 

Away from half filling, a nonzero chemical potential 
produces easy axis anisotropy due to the pseudospin de- 
pendence in the chemical potential term, thus favoring 
a CDW order. Furthermore, metallic physics of chargon 
excitations appears because Fermi-nesting is destroyed 
away from half filling. Then, the gapped phase of spinon 
excitations may be stable against gauge fluctuations ow- 
ing to the presence of gapless fermion excitations. [13] 
This anomalous metallic phase may be related with a 
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pseudogap phase in the context of high Tc superconduc- 
tivity. 

In tliis paper we developed a mean field theory tak- 
ing both electron excitations and order parameter fluc- 
tuations on an equal footing. The effective field theory 
consists of two parts. One is a fermion sector correspond- 
ing to the conventional Hartree-Fock or BCS mean field 
scheme. The other is a boson sector representing or- 
der parameter fluctuations that become fractionalized in 
the case of strong interactions, and are expressed as ro- 



tor variables. Softening of rotor fluctuations leads to a 
second order transition. Although our mean field analy- 
sis was performed in an insulating phase in the present 
paper, we believe that the present slave-rotor formula- 
tion would be more useful for studying quantum phase 
transitions in itinerant electrons because the presence of 
gapless fermion excitations can make the gapped phase 
of rotor excitations stable against gauge fluctuations. 

Extremely helpful discussions with Prof. J.-H. Han are 
appreciated. 
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